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Abstract 

Neutrino telescopes are moving steadily toward the goal of detecting astrophysi- 
cal neutrinos from the most powerful galactic and extragalactic sources. Here we 
describe analysis methods to search for high energy point-like neutrino sources us- 
ing detectors deep in the ice or sea. We simulate an ideal cubic kilometer detector 
based on real world performance of existing detectors such as AMANDA, IceCube, 
and ANTARES. An unbinned likelihood ratio method is applied, making use of the 
point spread function and energy distribution of simulated neutrino signal events 
to separate them from the background of atmospheric neutrinos produced by cos- 
mic ray showers. The unbinned point source analyses are shown to perform better 
than binned searches and, depending on the source spectral index, the use of energy 
information is shown to improve discovery potential by almost a factor of two. 



1 Introduction 

With the construction of IceCube at the South Pole [1] and of ANTARES in 
the Mediterranean Sea [2], together with existing R&D programs for a cubic 
kilometer array at these latitudes [3,4], neutrino astronomy is entering a very 
promising era. IceCube, when complete in 2011, will consist of up to 80 strings 
on a hexagonal grid with 124 m spacing, with each string holding 60 optical 
modules vertically spaced by 17 m. The instrumented part of the strings is 
deployed between 1450 and 2450 m deep in the ice. The experiment profits 
from experience acquired with the AMANDA detector, taking data since 1996 
and completed in 2000 with 19 strings containing 677 total optical modules 
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between 1500 and 2000 m below the ice surface. Studies on the optimal config- 
uration for a cubic kilometer array in the Mediterranean are ongoing, and the 
amount of photomultipliers (PMTs) generally considered is somewhat larger 
than for IceCube, around 6000 [5] and up to about 9000 [6] . 

The community is refining methods to detect low statistics signals amongst 
large backgrounds. The expected background from atmospheric neutrinos in 
a cubic kilometer detector is of the order of 50 000 upgoing events per year 
after selection criteria guaranteeing good angular resolution and rejection of 
misreconstructed cosmic ray muons. Point-like signals of few events need to 
be singled out among this large number of background events. Two features 
distinguish signal from the background: 

• The angular distribution. The signal would cluster around the direction of 
the neutrino source (assumed here to be point-like) with a spread depending 
on the detector angular resolution. Angular resolution is limited by detector 
geometry and by the propagation characteristics of light in the medium, 
specifically by photon scattering and absorption. The pointing accuracy for 
astrophysical sources is also limited by the kinematic angle between the 
parent neutrino and the muon. 

• The energy distribution. The differential energy spectrum of the signal ex- 
pected from Fermi acceleration mechanisms is close to E -2 , harder than that 
of atmospheric neutrinos due to the showering process in the atmosphere. 
The differential spectrum of atmospheric neutrinos approximately follows a 
power law of E~ 3 - 7 above 100 GeV. 

Other signatures may be used, including time dependencies of emissions mea- 
sured in other detectors such as correlations with gamma ray bursts or TeV 
gamma ray flares. However, we focus on steady emissions of neutrinos with 
time. The methods we have implemented exploit the two features listed above. 
We show that unbinned methods based on the likelihood ratio hypothesis test 
perform better than methods based on angular bins, and we show that the 
introduction of energy dependent information, e.g. the number of hit PMTs, 
helps discriminate signal and allows an energy spectrum reconstruction even 
when few events are detected on top of the background. Other unbinned meth- 
ods have been studied and developed by other authors [7,8,9,10]. 

Sec. 2 describes the simulation we use to generate realistic samples of signal 
and background events in a cubic kilometer detector. Sec. 3 describes the 
unbinned method based on a likelihood ratio analysis, comparing a signal 
plus background hypothesis to a background only hypothesis. This method 
has been applied to IceCube data for 9 strings and to 2005-6 data of the 
AMANDA-II detector [11,12]. In Sec. 4 we describe the performance of the 
methods and show results in terms of discovery potential. We also emphasize 
the importance of using energy related information to increase the discovery 
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Fig. 1. Simulated detector consisting of 81 strings, each containing 60 optical mod- 
ules. 

potential and show the ability to determine the spectral index of the neutrino 
source. The method is then compared to a search using angular bins, more 
traditionally applied in neutrino astronomy [13,14,15]. 



2 Simulation of a data sample of atmospheric neutrino background 
and point source signal 

We wish to compare several neutrino point source search methods and draw 
general conclusions on the discovery potential for cubic kilometer scale neu- 
trino telescopes under construction in the ice and under study in the sea water. 
We have performed a realistic simulation of both atmospheric neutrino events 
and signal events from an astrophysical neutrino point-like source. This is 
accomplished by a detailed detector volume simulation with simplifications 
expected to have negligible impact on the comparison. 

The simulated detector consists of 4860 optical modules arranged in 81 strings. 
The strings are evenly distributed on a square 9x9 grid with 125 meters 
separating nearest neighbors, and the modules on each string are vertically 
separated by 16.67 m, shown in Fig. 1. The modules are simulated as contain- 
ing a downward looking 10 inch photomultiplier with 20% quantum efficiency 
A multiplicity trigger requires that at minimum 14 modules register a pho- 



ton from an event. We neglect trigger time windows and photon hits from 
photomultiplier dark noise, since hits unrelated to the track can presumably 
be removed with coincidence requirements, filtering strategies, and topological 
cuts. We assume the detector is located at the South Pole. A different location 
would imply a different visibility of sources in the sky, since at the South Pole 
half of the sky is always visible while the other half is inaccessible. At other 
latitudes the detector is always blind to less than one half of the sky, always 
sensitive to a region of similar size in the opposite hemisphere, and sensitive 
to the remainder of the sky for a fraction of the day. 

We have prepared an algorithm to simulate the muon neutrino interactions 
in a volume and propagate the secondary muon in the ice, producing the 
light detected by the PMTs. For neutrino generation, we use an updated 
version of the simulation described in [16]. We use the more recent CTEQ6 [17] 
structure functions for the deep inelastic cross section of muon neutrinos and 
simulate the Earth density profile in [18] to account for the absorption of high 
energy neutrinos. The muon is propagated to the instrumented region using 
the MUM propagation code [19]. We model the muon energy loss within the 
detector as dE/dx = a + bE, where a = 0.268 GeV/m accounts for ionization 
energy losses, and b = 0.00047 m^ 1 for bremsstrahlung, pair production, and 
photonuclear interactions [20]. The number of PMTs recording light, or 'hit', 
increases with the energy loss rate, and thus the energy, of the muon. Since the 
energy dependent processes are stochastic, treating them as continuous may 
overestimate the number of modules hit by a small amount at PeV energies. 
However, high energy muons where this effect is significant are very bright and 
tend to produce photon hits in many more modules than the trigger threshold 
of 14, so we conclude that the impact on event triggering is negligible. 

Photon propagation in the detector medium has a significant impact on the 
detector response. We assume a homogeneous detector medium in which the 
photon density at the sensor can be described as a simple function with respect 
to the distance of the muon track. Photons are propagated with an effective 
scattering length of 21 m and an absorption length of 120 m, which are typical 
values for South Pole ice [21]. Fig. 2 shows the assumed photon density for 
a minimum ionizing muon as a function of distance for the case of ice. For a 
muon with arbitrary energy loss dE/dx, the photon density is scaled by the 
equivalent energy loss of a minimum ionizing muon, i.e. scaled by dE/dx / 
(0.268 GeV/m). The number of photoelectrons recorded from each PMT is a 
Poisson random variable with mean equal to the product of this photon den- 
sity, PMT photocathode area, and PMT quantum efficiency For each PMT, a 
hit is determined by randomly sampling this Poisson probability of observing 
at least one photoelectron. 

We have generated 10 11 upgoing neutrino events isotropically entering the 
Earth. The events are drawn from an E -14 energy spectrum between 10 and 
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Fig. 2. Photon density as a function of the radius from a minimum ionizing muon 
track in ice with an effective scattering length of 21 meters and absorption length 
of 120 meters. 

10 9 GeV. This hard spectrum allows efficient generation of high energy events. 
Events passing the trigger threshold of > 14 hit PMTs are kept, resulting in the 
neutrino effective area shown in Fig. 3. The neutrino effective area represents 
the equivalent detector area for a hypothetical instrument with 100% efficiency 
for detecting the passage of neutrinos. The effective area is much smaller 
than the dimensions of the detector due to the small cross section of neutrino 
interaction and, for larger zenith angles and high energies, neutrino absorption 
on transit through the Earth. It is a useful parameter for determining event 
rates and making comparisons between experiments. The event rate for a 
neutrino model predicting a flux IF d ^ n is given by 



Using the atmospheric neutrino flux of Barr et al. [22], we find approximately 
134 000 atmospheric neutrino events per year at trigger level with a median 
energy of 670 GeV and a slight density dependence on zenith angle, and thus 
declination, of ±15%. This number of events will be reduced however by topo- 
logical and reconstruction cuts necessary to eliminate the large background 
of misreconstructed downgoing muons from cosmic ray air showers. We esti- 
mate 50% of atmospheric events will be lost to reject this background; thus we 
choose 67 000 events from the atmospheric neutrino sample for our data set, 
roughly representing a year of data. Cuts applied to a real data sample may 
have some small energy dependence, i.e. efficiency may drop for low energy 
events, but we do not consider this effect. We generate astrophysical signal 
neutrinos similarly; however, they are simulated at a fixed declination and 
weighted differently than atmospheric neutrinos. We consider for the signal a 
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power law neutrino spectrum with spectral indices ranging from 1.5 to 4.0. 
The number of hit modules distribution for atmospheric neutrinos and several 
signal spectral indices is shown in Fig. 4. 



Angular reconstruction errors for a cubic kilometer detector in ice are shown to 
be ~0.7° for neutrinos of energies > 1 TeV [1]. We simulate this reconstruction 
error by adding to the muon direction a space angle error randomly sampled 
from a polar Gaussian distribution with standard deviation 0.7°. Considering 
both this reconstruction error and the neutrino-muon vertex angle, the result- 
ing median angular resolution for an E -2 signal is 0.86° and is larger for softer 
spectra. A comparison to a detector with angular resolution 0.2°, which may 
be achieved with a detector in sea water, is shown in Sec. 4. 



Muon energy reconstruction errors in neutrino telescopes are typically of the 
order of 0.3 in log 10 E v above a few TeV [23,24,25], limited by the stochastic 
nature of muon energy losses. To simulate this effect, we assign to each event 
a reconstructed error sampled from a Gaussian distribution of width 0.3 in 
the logarithm of energy. The reconstructed energy for atmospheric neutrino 
events and several power law spectra is shown in Fig. 4. Energy resolution 
degrades below a few TeV; thus the energy resolution we assign to such events 
is unrealistically accurate. However, the power to detect astrophysical sources 
resides in the ability to discriminate high energy neutrinos from lower energy 
atmospheric neutrino background, so we conclude overestimating the energy 
resolution of such lower energy events does not affect the result significantly. 
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Fig. 4. Number of hit modules distribution (left) and reconstructed muon energy 
distribution (right) for several neutrino spectra. 

3 The unbinned likelihood ratio method 

In the context of a muon neutrino point source search, the data from a neutrino 
telescope consists of a set of muon events spread over the sky, each with 
reconstructed direction (declination and right ascension), energy, and time. 
While event time is useful in searches for short neutrino bursts or periodic 
neutrino emission, we focus on searches for continuous neutrino emission and 
disregard event time. The vast majority of events are muons produced by 
atmospheric neutrinos. At any celestial direction, the data can be modeled by 
two hypotheses: 

• Hq: The data consists solely of background atmospheric neutrino events. 

• Hs'- The data consists of atmospheric neutrino events as well as astrophys- 
ical neutrino events produced by a source with some strength and energy 



The likelihood of obtaining the data given each hypothesis is calculable, and 
the ratio of likelihoods, or equivalently the log of the likelihood ratio, serves 
powerful test. We define our test statistic 



Larger values of A indicate the data is less compatible with the background hy- 
pothesis H . The probability density functions P(Data\H ) and P(Data\Hs) 
are calculated using knowledge of the spatial and energy distribution of back- 
ground and astrophysical neutrino events. 



spectrum. 




(2) 
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Suppose we wish to test the existence of a source at a known direction x s using 
a set of data events, each with reconstructed direction Xi and reconstructed 
energy E^. Events reconstructed outside a declination band centered at x s 
with a width several times the detector resolution are unlikely to have been 
produced by a source at x s and can be disregarded, leaving N events in the 
band. Each event inside the band is assigned a source probability density 
corresponding to the probability of the event belonging to a source at x s . The 
source is assumed to emit neutrinos according to an E~ 7 power law energy 
spectrum. The source probability density is the product of a spatial density 
function describing the potential of an event reconstructed with direction x*i 
to have true direction x s and the probability of observing reconstructed muon 
energy E i given source spectral index 7: 

Si{x h x s ,E hl ) = N{3i\£ a ) ■ J P{Ei\E v )P{E v \i)dE v . (3) 



The spatial probability density component N(xi\x s ) can be obtained directly 
from event reconstruction when maximum likelihood reconstruction techniques 
are used [26]. The structure of the reconstruction likelihood near the most 
likely direction Xi provides an event by event estimate of reconstruction un- 
certainty and yields a Gaussian spatial probability density profile: 

Mi{xi\x a ) = 7 r^e~ l ' El ^ L , (4) 



where I X £ X ^ I IS the space angle difference between source and reconstructed 
event directions, and a is the reconstruction error estimate. The vertex angle 
between the neutrino and muon is neglected in Mi{xi\x s ). This angle is negli- 
gible compared to reconstruction error for cubic kilometer neutrino telescopes 
in ice with resolution ~0.7° and high energy threshold of ~100 GeV. Addition 
of the vertex angle error is discussed in section 4 for cubic kilometer telescopes 
in sea water with resolution ~0.2°. The integral over E v is precomputed from 
detector Monte Carlo. The resulting tables of P(Ei\j) for 1.0 < 7 < 4.0, shown 
for several indices in Fig. 4, are computed in steps of 0.01 in 7, interpolated 
linearly, and stored for reference Q The resulting source probability density is 

Si{x h x s ,E hl ) = -Le-^^P^) (5) 



and has value unity when integrated over solid angle and E^. The slight decli- 
nation dependence of the atmospheric neutrino background can be neglected 

1 For maximum accuracy, one may wish to tabulate P(Ei\^) with respect to zenith 
angle as well; however, we omit this step and achieve good results. 
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over the width of the band. The background probability density again depends 
on event energy and is then 



P{Ei\(j) a t m ) 



(6) 



band 



The probability density P{Ei\4> atm ) is precalculated as above assuming the at- 
mospheric neutrino flux of Barr et al. [22] . We assume the background is purely 
atmospheric neutrinos. If the background contains another component, for ex- 
ample high energy muons from cosmic ray air showers, the energy component 
of the background density should be modified accordingly: P (E i \(f) a t m , + <t> Bkgd) ■ 
The source and background densities are combined, and the likelihood is eval- 
uated over all events in the band: 



where n s describes the number of signal events present in the band. 

The fraction of signal events n s as well as the source spectral index 7 are 
not known and must be determined by maximizing the likelihood C. This 
is done by minimizing the quantity —log(C) using the MIGRAD minimizer 
available in MINUIT [27] with respect to the unknown quantities n s and 7 
and obtaining the best value of each parameter, n s and 7. The minimization 
procedure also incorporates a penalty factoi(3]for 7 > 2.7 to better discrim- 
inate astrophysical sources with hard spectral indices from the atmospheric 
background with 7 ~3.7. The original hypotheses can be written in terms of 
C: P(Data\H ) = C(x s ,0) and P(Data\H s ) = n s , 7). The test statistic 



Ignoring the factor sign(h s ), the test statistic A is never negative since C(x s , 0) 
is contained in the range of C(x s , n s , 7), of which C(x s , n s , 7) is the maximum. 
A downward fluctuation of the background may occur at x s which would be 
fit best as a source with negative number of events and a negative value of n s . 
Such a downward fluctuation also would have a large value of A, so sign(n s ) is 
used to separate negative and positive excesses. A full sky search is a simple 
extension of this single point search method and can be accomplished by 

2 The extension of 7 is limited above 2.7 during the minimization by a Gaus- 
sian likelihood penalty with a=0.2 in spectral index. This improves the discovery 
potential for hard astrophysical source spectra while remaining comparable to an 
unbinned method without energy for source spectra as soft as ~3.2 (see Fig. 8). 




(7) 



is 




(8) 
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performing the search on a grid of locations covering the sky. Finally, while it 
is preferable to use an energy estimation to maximize the power to discriminate 
astrophysical neutrinos from the background, it is possible to do the search 
without this information. The energy dependent terms are removed (i.e., set 
to one) from the signal and background probability densities, 

1 _\x i -x s \ 2 

Si(xi,x s ) = - — (9) 
Bi = ^— (10) 

1 ^band 



resulting in simpler expressions for L. The quantity —log(C) is minimized with 
respect to n s , and P(Data\H ) = £(x s ,0) and P(Data\H s ) = C(x s ,n s ). 



4 Results 



We apply the likelihood ratio method to data consisting of the 67 000 back- 
ground events described in section 2 and an added source at declination 48°. 
A grid of simulated source strengths and spectral indices is used, with source 
strength spanning - 100 signal events added to the sample and spectral index 
7 spanning 1.0 - 3.9 in increments of 0.1, resulting in a total of 3000 simulated 
combinations of source strength and spectral index. For each combination, 10 
000 experiments are done and the value of A is recorded for each. 10 7 trials 
are performed with background alone to evaluate the significance of observed 
values of A. The method is compared against a binned search with a circular 
bin centered at the source location. The bin radius is optimized by minimizing 
the number of E" 2 signal events necessary to achieve 5a significance in 90% of 
experiments. The optimal radius is found to be 1.35°, with a signal efficiency 
of 80% and background expectation of ~16 events/bin. Also, the method is 
compared against itself using only spatial information and neglecting event 
energy. An identical grid of simulated source strengths and spectral indices is 
used for both the binned method and likelihood method without energy. 

Fig. 5 shows a significance sky map with an added source at declination 5=48° 
and right ascension a=12h producing 15 events according to an E~ 2 energy 
spectrum. Fig. 6 illustrates the procedure used to determine significance and 
discovery potential. The integral distribution of A for background alone is 
produced at declination 5=48° and the values of A corresponding to 3<r (2.7 x 
10~ 3 ) and 5<T (5.7 x 10~ 7 ) integral probability are calculated. Fig. 6 also shows 
distributions of lambda with 8, 16, and 24 signal events added to background. 
Discovery potential at 5<r is the fraction of experiments with A exceeding the 
5cr threshold. Discovery potential is then computed in this fashion for each 
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Fig. 5. Significance sky map of 67 000 background events with an added source 
of 15 events distributed according to E~ 2 at declination 5=48° and right ascension 
a=12h 
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Fig. 6. Left: Integral A probability distribution for background at declination 5=48°. 
Shaded regions show statistical uncertainty. The 2 horizontal lines indicate the val- 
ues of the probabilities corresponding to 3a and 5a confidence levels. Right: Dis- 
tribution of A for background at declination 5=48° and with 8, 16, and 24 added 
signal events distributed according to an E -2 energy spectrum. 

number of signal events and each spectral index at 5a confidence level. The 
number of detected events produced by a source of a given strength are Poisson 
distributed around a mean strength related directly to the source flux. The 
detection probability for a given mean strength is calculated by summing over 
the detection probabilities for all - 100 signal events and weighting by the 
appropriate Poisson probability. The detection probability at 5a confidence 
level as a function of mean source strength for an E~ 2 signal flux is shown 
in Fig. 7. As can be seen, the likelihood method requires approximately half 
the signal flux needed by the binned method to reach a similar detection 
probability. To reach the same detection probability without using energy, the 
likelihood search still requires ~10% less signal flux than the binned method. 
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Fig. 7. 5cr confidence level detection probability vs. Poisson mean number of source 
emitted events added to the background sample of 67 000 events, representing one 
year of data from a km 3 scale neutrino telescope. The solid line indicates discovery 
potential for this method, the dashed line for this method without using energy 
information, and the dotted line for the binned method. The events follow an E~ 2 
energy spectrum. 



In Fig. 8 we show the 50% detection probability at 5a confidence level as a 
function of source spectral index for the three methods along with the values 
obtained from a similar analysis with a detector capable of a muon angular 
resolution of 0.2°, similar to the resolution possibly achieved by km 3 detectors 
in water [2,5]. With an angular resolution 0.2°, the neutrino-muon vertex angle 
can no longer be neglected in the likelihood method. The resolution of 0.2° is 
convoluted with the neutrino-muon vertex angle as a function of reconstructed 
muon energy for an E~ 2 neutrino spectrum in bins of 0.2 in logioE. The con- 
voluted resolution is used to compute Ni(xi\x s ) in the likelihood method. The 
bin size used by the binned method in water is reoptimized and found to be 
0.45°. Both detectors, ice and water, are simulated at the South Pole to avoid 
differences in visibility as a function of declination. 

For a source with an E™ 2 energy spectrum, the binned method requires 30 
signal events on top of background for a 50% chance of 5a detection, while 
the likelihood method requires 26 events without using energy or only 14 
events with energy. For all methods, harder spectral indices require fewer 
events for 5a detection. This trend is caused by two factors: Higher energy 
events have a smaller neutrino-muon vertex angle and thus better resolution, 
and high energy astrophysical neutrinos produce muons at the detector with 
much higher energy than typical atmospheric neutrinos. The binned method is 
affected only by changes in the vertex angle, while the likelihood method using 
event energy is affected by both factors and has a stronger response to spectral 
index. In ice, the method improves from only slightly better than binned 
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Fig. 8. Poisson mean number of signal events required for 50% detection probability 
at 5a confidence level vs. source spectral index for the unbinned method described 
in the paper using energy information in water (solid lower line) and in ice (solid 
upper line) , not using energy information in ice (dashed line) and for a binned search 
for a detector in water (dotted lower line) and in ice (dotted upper line). For the 
detector in ice we assumed an angular resolution of 0.7°, and we assume a resolution 
of 0.2° in water. For both detectors, a background of 67 000 events is used. 

methods for soft spectral indices greater than ~3 to more than a factor of two 
improvement for hard spectral indices less than ~2. A possible way to improve 
the binned method is application of a cut on muon energy, keeping high energy 
events and reducing background, but such a cut must be optimized assuming a 
specific source spectral index and is generally undesirable. Optimizing search 
bin size for a specific spectral index is similarly undesirable. The likelihood 
method requires no a priori assumption of signal spectral index to utilize event 
energy in discriminating signal from background. 

Assuming a similar atmospheric neutrino background, the number of events 
required for 5cr discovery in sea water is a factor of ~2 less than for ice if 
a resolution of 0.2° is achieved. For a detector in the Mediterranean, more 
than a year of livetime may be required to reach this background rate for 
declinations with less than 100% exposure. With high statistics, significance 
improves inversely with resolution, so one might expect a detector in ice with 
a resolution of 0.7° to require a factor of 3.5 more events than one in sea water 
with a resolution of 0.2° to achieve a similar significance. However, in this 
analysis the factor is smaller because event energy discriminates signal from 
background, and background rates are relatively low. 

Finally, since source spectral index is a free parameter in the method and fit- 
ted to the most likely value, the method provides an estimate of the spectral 
index. Fig. 9 summarizes this capability for source intensities of 15 events and 
50 events following an E~ 2 energy spectrum at declination 5=48°. The distri- 
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Fig. 9. Confidence boundaries are shown in the spectral index vs. source strength 
parameter space for 15 signal events (left) and 50 signal events (right) for 67% and 
90% confidence levels. 

bution of —2 • log[C(n s , j)/C(n s , 7)] is approximately given by a chi square 
with two degrees of freedom. Confidence contours are produced by comparing 
—2 • log likelihood ratio values for points on the plane of source strength and 
spectral index to chi square values corresponding to 67% and 90% confidence 
level. The coordinates of the best fit (the largest value of A) and true point are 
indicated. Naturally, the confidence regions shrink for larger signal strengths, 
indicating better estimation of both parameters. Good spectral index recon- 
struction is achieved with few events from a source despite limited energy 
resolution on an event by event basis. This is possible because neutrino tele- 
scopes are capable of detecting neutrinos over many energy decades, offering 
a large lever arm to discriminate energy spectra. 

For a detector at the South Pole, background event density and angular res- 
olution do not strongly depend on declination, so the results shown at 5=48° 
generalize for declinations above 5=0°. A detector in the Mediterranean will 
have stronger declination dependence of background event density due to ex- 
posure differences, so the number of events required for 5<r detection becomes 
smaller for declinations with reduced exposure. 



5 Conclusions 



We have shown that likelihood methods improve discovery potential using dif- 
ferences in both the angular distribution of events and energy spectra between 
the background and signal hypotheses. Furthermore, it has been shown that 
10% less flux is required compared to methods using predefined angular search 
bins. The improvement is more significant when an energy related variable is 
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used to distinguish hard astrophysical neutrino spectra from the softer spec- 
trum of atmospheric neutrinos. In this case, for E~ 2 energy spectra only about 
half of the flux is required to achieve 5a compared to binned methods. 



References 



[1] J. Ahrens et al, The IceCube Collaboration, Astropart. Phys. 20 (2004) 507; 
http://icecube.wisc.edu, 

[2] The ANTARES Collaboration, A Deep Sea Telescope for High Energy Neutrinos 



- Proposal for a 0.1km2 detector by the ANTARES collaboration, eprint: astro- 
ph/9907432||http://antares.in2p3.fr 



[3] The KM3NeT Collaboration, Design Study for a Deep Sea Facility in the 
Mediterranean for Neutrino Astronomy and Associated Sciences, Proposal, 
available at |http: / /www.km3net.org/proposal/KM3NeT-FINAL-PROPOSAL^ 
ShorterVersionForWebsite.pdf and http:/ /km3net.org. 



[4] |http: / /nemo.lns.infn.it 



[5] S. Aiello et al, Astropart. Phys. 28 (2007) 9. 

[6] J. Carr, D. Dornic, F. Jouvenot, G. Maurin for the KM3NeT Collaboration, 
Configuration studies for a cubic-kilometer deep-sea neutrino telescope 
- KMSNeT - with NESSY, a fast and flexible approach, proc. of 
Int. Cosmic Ray Conf. (ICRC07), Merida, July 2007, available at: 



http://www.km3net.org/talks/2007/icrc2007/icrc0859.pdf. 



[7] T. Neunhoffer and L. Kopke, Nucl. Instr. Meth. A 558 (2006) 561. 

[8] A. Heijboer, Track reconstruction and pointsource searches with Antares, PhD 
Thesis (2004), University of Amsterdam, available at 
http: / /antares. in2p3.fr/Publications/index.html#thesis , 

[9] C. Finley and S. Westerhoff, Astropart. Phys. 21 (2004) 359. 

[10] J. A. Aguilar Sanchez for the ANTARES Collaboration, Expected discovery 
potential and sensitivity of the ANTARES neutrino telescope to neutrino point- 
like sources, proc. of 30* h Int. Cosmic Ray conference (ICRC2007), July 2007, 
Merida, Mexico and e-Print: [arXiv:0710.0 252 [astro-ph]. 

[11] C. Finley for the IceCube Collaboration, Nine- string IceCube point source 
analysis, proc. of 30 t/l Int. Cosmic Ray conference (ICRC2007), July 2007, 
Merida, Mexico and e-Print: |a7Xiv:071 1.0353] [astro-ph], pp. 107-110. 

[12] J. Braun for the IceCube Collaboration, Neutrino point source search strategies 
for AMANDA-II and results from 2005, proc. of 30 t/l Int. Cosmic Ray conference 
(ICRC2007), July 2007, Merida, Mexico and e-Print: |arXiv:0711.0353| [astro- 
ph], pp. 99-102. 



15 



[is; 

[14 
[15 

[16 

[i7; 
[is; 

[19 
[20 

[21 

[22 
[23; 

[24 



[25 



[26; 
[27 



M. Ambrosio et al, Astrophys. J. 546 (2001) 1038. 
A. Achterberg et al, Phys. Rev. D75 (2007) 102001. 

K. Abe et al High energy neutrino astronomy using upward-going muons in 
Super-Kamiokande-I, Astrophys. J. 652 (2006) 198. 

E. Bugaev, T. Montaruli, Y. Shlepin and I. Sokalski, Astropart. Phys. 21 (2004) 
491. 

J. Pumplin et al, JHEP02 (2006) 032, structure functions available at 
|http://user .pa.msu.edu/wkt/cteq/cteq6/cteq6pdf.htmrj 

A. Dziewonski, Encyclopedia of Solid Earth geophysics, ed. by D. E. James 
(Van Nostrand Reinhold, New York, 1989), p. 331; we used the formula from 
this work also cited in R. Gandhi, C. Quigg, M. Reno and I. Sarcevic. 

I. Sokalski, E. Bugaev and S. Klimushin, Phys. Rev. D64 (2001) 074015. 

D. Chirkin and W. Rhode, Muon Monte Carlo: a high-precision tool for muon 
propagation through matter, e-Print:hep-ph/0407075. 

M. Ackermann et ah, Optical Properties of deep glacial ice at the South Pole, 
J. Geophys. Res. Ill (2006) D13203. 

G. Barr et al, Phys. Rev. D70 (2004) 023006. 

J. Ahrens et al, the AMANDA Collaboration, Nucl. Instrum. & Meth. A524 
(2004) 169. 

J. D. Zornoza and D. Chirkin for the IceCube Collaboration, Muon energy 
reconstruction and atmospheric neutrino spectrum unfolding with the IceCube 
detector, proc. of 30 t/l Int. Cosmic Ray conference (ICRC2007), July 2007, 



Merida, Mexico and e-Print: arXiv:0710.0252 [astro-ph]. 



A. Romeyer et al. for the ANTARES Collaboration, Muon energy reconstruction 
in ANTARES and its application to the diffuse neutrino flux, proc. of 28 tft 



Int. Cosmic Ray Conf. (ICRC2003), Tsukuba, Japan, July 2003, e-Print: astro- 
|ph/0308074| 

T. Neunhoffer, Astropart. Phys. 25 (2006) 220. 

F. James, MINUIT: Function Minimization and Error Analysis Reference 



Manual, available at http://wwwasdoc.web.cern.ch/wwwasdoc/mimit7J 



16 



